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Abstract 

By modeling a linear, anisotropic and inhomogeneous magnetodi- 
electric medium with two independent set of harmonic oscillators, elec- 
tromagnetic field is quantized in such a medium. The electric and mag- 
netic polarizations of the medium are expressed as linear combinations 
of the ladder operators describing the magnetodielectric medium. The 
Maxwell and the constitutive equations of the medium are obtained as 
the Heisenberg equations of the total system. The electric and mag- 
netic susceptibilities of the medium are obtained in terms of the ten- 
sors coupling the medium with the electromagnetic field. The explicit 
forms of the electromagnetic field operators are obtained in terms of 
the ladder operators of the medium. 
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1 Introduction 



One method to quantize the electromagnetic field in the presence of an ab- 
sorptive medium is known as the Green function method PQ-[E]. In this 
method by adding the noise electric and magnetic polarization densities to 
the classical constitutive equations of the medium, these equations are taken 
as definitions of electric and magnetic polarization operators. The noise 
polarizations are related to two independent sets of bosonic operators. Com- 
bination of Maxwell and the constitutive equations in frequency domain, 
gives the electromagnetic field operators in terms of the noise polarizations 
and classical Green tensor. The commutation relations are imposed on the 
bosonic operators such that the commutation relations between electromag- 
netic field operators in the magnetodielectric medium become identical with 
those in free space. Another interesting quantization scheme of electromag- 
netic field in the presence of an absorptive dielectric medium is known as 
the damped polarization model which is based on the Hopfield model of 
a dielectric [DJ, where the polarization of the dielectric is represented by a 
damped quantum field [ID] . In the damped polarization model [H]-[l3] , the 
electric polarization of the medium is represented by a quantum field and 
absorptivity character of the medium is described by interaction between the 
polarization with a heat bath containing a continua of harmonic oscillators. 
In this method a canonical quantization is formulated for electromagnetic 
field and the medium. The dielectric function of the medium is obtained 
in terms of coupling function of the heat bath and the electric polarization, 
such that it satisfies the Kramers- Kronig relations [H] . Recently Raabe and 
et al [T5] have represented a unified method to quantize electromagnetic field 
in the presence of an arbitrary linear medium based on a general nonlocal 
conductivity tensor and using a single set of appropriate bosonic operators. 
This formalism recovers and generalizes the previous quantization schemes 
for diverse classes of linear media. In particular the quantization of electro- 
magnetic field in the presence of a magnetodielectric medium is a limiting 
case for weakly spatially dispersive medium in this scheme. In the present 
work, we generalize our approach [TBI [T7] to quantizing electromagnetic field 
in an anisotropic magnetodielectric medium with spatially and temporarily 
dispersive property. In this case, the electric and magnetic polarizations are 
dependent on the macroscopic electric and magnetic fields inside the medium 
in a non local way with respect to both position and time. The electric and 
magnetic polarization densities of the medium are defined as linear combi- 
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nations of the ladder operators of the medium. The coefficients of these 
linear expansions are coupling tensors which couple the medium with elec- 
tromagnetic field. The electric and magnetic susceptibility tensors of the 
medium are obtained in terms of the coupling tensors. By using a Hamilto- 
nian in which the electric and magnetic polarizations minimally couple to the 
displacement and magnetic fields respectively, both Maxwell and the consti- 
tutive equations of the medium can be obtained as the Heisenberg equations 
of the total system. Finally, using the Laplace and Fourier transformations, 
we obtain the space-time dependence of electromagnetic field operators in 
terms of the annihilation and creation operators of the oscillators modeling 
the medium. 



2 Quantization scheme 

In order to quantize electromagnetic field in the presence of an anisotropic 
magnetodielectric medium, we enter the medium directly in the process of 
quantization modeling it by two independent set of harmonic oscillators 
which we call them E and M quantum fields [16] - The E and M fields de- 
scribe polarizability and magnetizability of the medium, respectively. This 
means that in our approach the electric and magnetic polarization densities 
of the medium are defined respectively as linear combinations of the ladder 
operators of the E and M quantum fields. We use the Coloumb gauge in 
this quantization scheme and do the quantization in unbounded space and 
in the absence of external charges. Generalization of the quantization inside 
a cavity with a definite volume and with a known boundary conditions and 
or in the presence of external charges is straightforwardly [17] . 

Applying the Coloumb gauge, the quantum vector potential can be ex- 
panded in terms of plane waves as 



A(f,t) = Jd 3 kj2 



{\ 2(2tt)% u; 



'' «uW^+«Lw e_ 1^ (i) 



fcA v 

k 



where uj^ = c\k\, Eq is the permittivity of the vacuum and eg A , (A = 1,2) 
are unit polarization vectors which satisfy the orthogonality relations 

& kX -k = 0, (2) 
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Operators a% x (t) and al (t) are annihilation and creation operators of the 
electromagnetic field and satisfy the following equal time commutation rules 

Mt) ia ]f x Xt)) = 6(%-9)6 w . (3) 

Quantization in coloumb gauge usually needs resolution of a vector field in its 
transverse and longitudinal parts. Any vector field F(f) can be resolved in 
two components, transverse and longitudinal components which are denoted 
respectively by F 1 - and F« pT8| [19] . The transverse part satisfies the coloumb 
condition V ■ F = and the longitudinal component is a conservative field 
V x F" = 0. In the absence of external charges the displacement field is 
purely transverse and can be expanded in terms of the plane waves as 

D(r, t) = - lEo J d 3 qJ2 ~ «SX(*)^*£X- ( 4 ) 

where /io is the permeability of the vacuum. From ([3]) the commutation rela- 
tions between the components of the vector potential and the displacement 
field clearly are 

[Ai(r,t), -Dj(?,t)\ = ih5tj(r- p), (5) 

where 5y(r — r') = J d z ke lk '^ r ~ r '^ ] (8ij — is the transverse delta func- 
tion. 

Now we enter the medium in the process of quantization taking its con- 
tribution in Hamiltonian of the total system ( electromagnetic field plus 
medium) as the sum of the Hamiltonians of the E and M quantum fields 

Hd = H e + H m , 

H e(t) = I d3 Q I d 3 k huj^ di{k,q,t)d u {k,q,t), 

3 

H m (t) = ^ f d 3 qj d 3 k hu,:bt(k,q,t)K(k,q,t). (6) 

where H e and H m are the Hamiltonians the E and M fields respectively. 
Quantum dynamics of a dissipative harmonic oscillator interacting with an 
absorptive environment can be investigated by modeling the environment by 
a continuum of harmonic oscillators [20]- [28]. In the case of quantization of 
electromagnetic field in the presence of a magnetodielectric medium, electro- 
magnetic field is the main dissipative system and the medium plays the role 
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of the absorptive environment. Here electromagnetic field contains a contin- 
uous set of harmonic oscillators labeled by k and v — 1,2. Therefore to each 
harmonic oscillator of the electromagnetic field a continuum of oscillators 
should be corresponded. In the present scheme, to each harmonic oscilla- 
tor of the electromagnetic field labeled by k and is, we have corresponded 
two continuous set of harmonic oscillators defined by the ladder operators 
d u (k,q,t), dl(k,q,t) and b u (k,q,t), bl(k,q,t) which are to describe the elec- 
tric and magnetic properties of the medium respectively and satisfy the equal 
time commutation relations 

[du(k, q, t),dl(k', q', t)] = 5 w ,8{k - k')5(q- <?), 
[b v (k,q,t),bl(9,q\t)]^5 wl S(k-k')S(q-q"). (7) 

Summation on v = 3 in ([6]) is necessary, because as we will see the polariza- 
tion densities of the medium are defined in terms of the ladder operators of 
the E and M fields and in opposite to the vector potential the polarization 
densities are not purely transverse. In relation (Q ujj: is called dispersion rela- 
tion of the magnetodielectric medium and can be chosen simply as uj^ = c\k\ 
[T6l [T7| . It is remarkable that, although the medium is anisotropic in its 
electric and magnetic properties, we do not need take the dispersion relation 
as a tensor. In fact the dispersion relation has not any physical meaning 
here, and the Hamiltonian of the medium as is merely a mathematical 
institution so that gives us the correct form of the equation of motion of the 
total system, that is the Maxwell and the constitutive equations of the mag- 
netodielectric medium. The anisotropic behavior of the medium is expressed 
merely as definitions of the electric polarization density, denoted by P and 
magnetic polarization density, denoted M, in terms of the ladder operators 
of the E and M fields as follows 



5 < fq_ 

x [d„(k,q,t)e^ p + H.c.}vi(q), (8) 



Pi(r,t) = E / -T== I [ d 3 r%(^,r,P) 



Mi(f,t) 



•E 



d 3 q 



d s k / cPr'gijfafr^r') 



'(27T) 3 

x [b u (k,q,t)e l * p -H.c]si(q), 



(9) 
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where 



u=l,2 



Su{q) = q x e u $, is = 1,2 



Mo) = M$) = q 



\q\ 



(10) 
(11) 

(12) 



In (|Sj) and the fij(uj^ } f, r') and gij(u)%, r, r'), are some real valued tensors 
and are called the coupling tensors of electromagnetic field and the medium. 
The coupling tensors play the key role in this method and are a measure for 
the strength of the polarizability and magnetizability of the medium macro- 
scopically, so that we will see that the electric and magnetic susceptibility 
tensors of the magnetodielectric medium are obtained in terms of these cou- 
pling tensors. Also, the explicit forms of the noise polarization densities are 
obtained in terms of the coupling tensors and the ladder operators of the E 
and M fields at t = 0. It can be shown that when the medium tends a non 
absorbing one, the noise densities tend to zero and this quantization scheme 
reduce to the usual quantization in these media |16j . 

Now the total Hamiltonian, i.e., electromagnetic field plus the E and M 
quantum fields can be proposed as one of the following forms 



H(t) 



d 3 r 



[D-Pf [Vxl- ft ilf 



2e 



2^c 



+ H e + H ri 



(13) 



H(t) 



d A r 



[D - Pf [V x Af 



2a*c 



V x A-M 



(14) 



in which the electric and magnetic polarizations interact minimally with the 
displacement field and the magnetic field, respectively. Using both of the 
Hamiltonians ffl3|) and (|T4"|) give us correctly the Maxwell and constitutive 
equations of the manetodielectric medium as the Heisenberg equations of the 
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total system. Here we use the Hamiltonian (TT5]) since it is easier to solve the 
coupled Maxwell and constitutive equations especially when the medium is 
translation invariant in its electric and magnetic properties. 



3 Heisenberg equations 

3.1 Maxwell equations 

The Heisenberg equations for the fields A and D are 

^m = ^ m= . m^m t (15) 

= ' { H, n { f, t)] = VxV ;/>-'> - V x M(f, ( ), (16) 

where P 1 - is the transverse component of P. If we define the transverse 
electric field E- 1 , magnetic induction B and magnetic field H as 

E ± = B = VxA, H= — -M. (17) 

at 

The Heisenberg equations (|T5l) and (fl6l) can be rewritten as 

D = e E ± + P ± , (18) 
3D 

^^Vx«, (19) 

which are respectively the definition of the displacement field and the macro- 
scopic Maxwell equation in the absence of external charges. In the absence of 
external charges. In coloumb gauge we can take the longitudinal component 
of the electric field as iP = — According to the definitions ([T7]) we have 

an 

Vx£=--. (20) 
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3.2 Constitutive equations of the magnetodielectric medium 

Using commutation relations ([7j) we easily find the Heisenberg equations for 
operators d u (k,q,t) and b u (k,q,t) as 

d v (k, q, t) = ^[H, d u (k, q, t)} = -uo^ d u (k, q, t) 
+ ,\ = I d 3 r' J d 3 r" e^'fijiwt,?,-?') E l (r',t) 

(21) 



b u (k,q,t) = - [H,b u (k,q,t)] = -uv% b„(k,q,t) 

+ 1"° I d 3 r' f d 3 r" e-^>£,rV) H\r>,t) si(q). 
hJ(27r) 3j J 



(22) 



These equations can be solved formally as 



dJk, q, t) = dJk, q, 0)e-^E* + — , / dt' e^S^ 

x / d 3 r' I d 3 r" e^f ij (tJ S ,?,r*)E i (?J)vi@, 



(23) 



Uk, q, t) = bJk, q, 0)e-^* + /* dt' e -^(*-*') 

nJ(2vr) 3 -Jo 



x y rfVy dV' e" 4 ^ 9lJ {uJ h r'^')H\r',t')si{q). 

(24) 

Substituting of (123!) in (jSJ) and (!24j) in (jUJ) and using the completeness relations 

3 3 

i/=i i/=i 

give us the constitutive equations of the magnetodielectric medium which 
relate the electric and magnetic polarization densities of the medium to the 
macroscopic electric and magnetic fields respectively as 

P,(f, t) = P Nl (f, t) + e jf " dt' J d 3 r' 4-(f, r>, \t\ - ±t'), (26) 
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Mi(r, t) = M m (f, t) + / ' dt' I d 3 r' X ™(r, r' , \t\ - t')W{r', ±f ), (27) 

where the upper (lower) sign corresponds to t > (t < 0) and E — — ^ — 
is the total electric field. The memory tensors 

x e (f,r',t) = 



8tc 
hc 3 e Jo 



duuj z sin u)t I d 3 r" 



f{u, f , r 7 ') • f\u , P , r 7 ') 



t > 

t < 
(28 



and 



X m (f,r',t) 



' 871-^0 f°° 



dijju) sino>t / d r 



g(o>, r , r") • g l {uj ,r' , r") 



t > 



£ < 



(29) 



are respectively the electric and magnetic susceptibility tensors of the mag- 
netodielectric medium which have been obtained in terms of the coupling 
tensors / , g and their transpositions /* , g l . For a medium with a definite 
pair of the tensors \ e an d X m it is possible to solve equations f[2"8"|) and fT2"9l 
in terms of the coupling tensors / and g using a type of eigenvalue problem 
[T5] . The operators Pn and in (126|) and (I2?|) are the noise electric and 
magnetic polarization densities and their explicit forms are given by 



P m (f,t) = £ / -4= I f d 3 r' fiA^fy) 



x 



(2tt) 3 

[d„(k,q,0)e- l ^ t+l * p + H.c] vi(q), 



(30) 
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M Ni {r,t) 



3 



d?q 



x 



[b u (k,q,0)e-^ t+ ^ p - H.c] 4(g). 



(31) 



From (1281) and (1291) it is clear that for a given pair of the susceptibility tensors 
X e and x" 1 , the solutions of the relations (l28]) and (129]) for the coupling tensors 
/ and g are not unique. In fact for a given pair of x e an d X m \ if the tensors 
/ and g satisfy equations (1251) and (1291) . then the coupling tensors 



(32) 



f\u,r,r') = J d 3 sf(u,r,s) ■ A\u,s,r') 
g'(u,r,r') — J d 3 s g(oo,f,s) ■ A*(o;,s, r') 

for any tensor A where satisfy the orthogonality relation 

d 3 r" A(u, a, r 7/ ) • A*(u, P, r 7 ') = J<5(s - P) 



(33) 



also satisfy equations (l2~g|) and (12^1) . Although this affect the space-time de- 
pendence of the noise polarizations and therefore the space-time dependence 
of the electromagnetic field operators but all of these are equivalent. This 
means that various choices of the coupling tensors / and g satisfying (1281) 
and (129|) . for a given pair of the susceptibility tensors \ e an d X n ' \ do not 
affect the commutation relations between the electromagnetic field opera- 
tors and hence the physical observables. This subject becomes more clear if 
we compute the commutation relations between the components of the tem- 
poral Fourier transformations of the noise polarizations Pn and using 
commutation relations (JTj) and obtain 



[P Ni {r,u) , Pj fj (?,u')] = 
[M m {f,u) , M j Nj (r',u')} 



he 



TV 



Im 



Xtj{r,r',u) 8(u-u/) 



_h_ 

/i 7T 



-Im 



X™(r,P,u)}5(u-uj') (34) 



where x e an d X m are respectively the temporal fourier transformations of 
the tensors x e an d X m ''■ The commutation relations (1341) are generalization 
of those in reference [8J for an anisotropic magnetodielectric medium with 
spatial-temporal dispersion. For a given pair of the tensors x e an d X m ; var_ 
ious choices of the coupling tensors / and g satisfying the relations (l2"gj) and 
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( 129)) . do not affect the commutation relations (I34p and therefore the commu- 
tation relations between the electromagnetic field operators. Therefore all of 
the field operators which are obtained with a definite pair of the susceptibil- 
ity tensors x e an d X m but with different coupling tensors satisfying (1281 and 

( 129) ) . are equivalent. 

It is clear from (I3"U|) and (|3"TT) that explicit forms of the noise polarization 
densities are known. Also, because the coupling functions /, g are com- 
mon factors in the noise densities Pn, Mn and the susceptibility tensors % e , 
X m , it is clear that the strength of the noise fields are dependent on the 
strength of x e an d X m which describe dissipative character of the magnetodi- 
electric medium. When the medium tends to a nonabsorbing one the noise 
polarization tends to zero and this quantization method is reduced to the 
quantization in the presence of a nonabsorbing medium [T6] . 

It should be noted that the time derivative of the polarization fields ^ 

and are continuous at time t = although the absolute value \t\ is 
appeared in the constitutive equations (126")) and (127)) . In fact an other solution 
for the Heisenberg equation (121)) can be written as 

dv(k, q, t) = 4 n (jfe, q)e^ + / C dt' e'^^ 

x J d 3 r' J We-Vfoiufr P, t')vi(q) (35) 

where d™(k, q) are some time independent annihilation operators which sat- 
isfy the same commutation relations (j7j). If we find d u (k, q, 0) from the equa- 
tion (135)) and substitute it in the noise polarization field P/v(r, t) given by 

( 130) ) . we deduce 

P Ni (r, t) = P™(r, t) + e f dt 1 f d 3 r%(f, f>, t - t')&(r>, t') 

J —OO J 

-so J%t' J d 3 r%(r,f',\t\-t')E j (r',±t'), (36) 

where the susceptibility tensor x e is given by (1281) and P™ is the same as 
Pjsf with the exception that the annihilation operators d u (k,q,0) should be 
replaced by c^ n (A;, q). Now sustituting P^ from (136)) into (1261) give us the 
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expression 

Pi{r, t) = P™(r, t) + e f dt' I dh'xW, ?, t - t')E\r>, f) 

J —oo J 

(37) 

for the polarization field P which is valid for both positive and negative times. 
From (1371) it is clear that ^ and accordingly the electromagnetic field oper- 
ators are continuous at t — 0. One can apply the constitutive equation (j3"7j) 
and similar equation for the magnetic polarization M and use the tempo- 
ral Fourier transformation, or apply the constitutive equations (1261) and (1271) 
and use the Laplace transformation, to solve the coupled constitutive and 
Maxwell equations. Here we prefer the later, since in this way it is easier to 
show the limiting cases in the absence of any medium or in the presence of 
a nonabsorbing medium [T6l [T7] . Using the forward and backward Laplace 
transformation and applying the constitutive equations fl26|) and (1271) , we can 
obtain the explicit forms of the electromagnetic operators for both negative 
and positive times which are continuous at t = 0[12j. 

4 Solution of the Heisenberg equation 

In this section we solve the coupled Maxwell and constitutive equations 
(|T8|) . (jT9l) . (|20|) . (1261) and (1271) for a translationally invariant medium, that 
is for a medium that its electric and magnetic susceptibility tensors are de- 
pendent on the difference f — r' . For such a medium we can obtain the 
explicit forms of the electromagnetic field operators easily for both positive 
and negative times by the spatial fourier transformation and the temporal 
Laplace transformation. In this section we solve the coupled Maxwell and 
constitutive equations for positive times using the forward Laplace transfor- 
mation. The solution of the Heisenberg equations for negative times can be 
found similarly from the backward Laplace transformation [12] . Let us define 

F_ for any vector field F(f, t) by 

F(k,p) = J d 3 r ^ dt F(r,t) e - lk ^ pt . (38) 

Now applying such a transformation on the both sides of the equations the 
constitutive equations ( fT8|) . ( 1261) and ( |27l) . together with the Maxwell equa- 
tions ( |T9l) and ( l20l) . and then their combination for a translationally invariant 
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medium we find 



— ik x E = -pp(k,p)H_ - n pM_ N (k, p) + R(k, 0), 
-ikx H = pi{k, p)E + pP N (k, p) - D(k, 0) , 



(39) 



where B_(k, 0) , D_(k, 0) are respectively the spatial fourier transformations of 
B(r, 0) and D(f, 0) and e = Eo(l+x e ) , p = po(l+X m ), are the transformation 
introduced in ( |38l) for permittivity and permeability tensors of the medium. 
The equations ( 139|) can be written in a compact form using a matrix notation 
as follows 



A(k,p) 



' E ' 




. K _ 





-pE N + D(k,0) 



(40) 



where A(k, p) is a 6 x 6 matrix defined by 

0(k) -pp(k, p) 



A(k,p) 



and 



0(k) 



pi(k,p) o(k) 

— tk 3 iki 



(41) 



ik 3 
— iki iki 



-ik\ 




(42) 



Finally, substituting P_ N (k,p), M N (k, p), D(k, 0) and B(k,0) in the right 
hand side of ( l4"0l) using the expansions (JT|) , (SI) , (1301 and ( |3T1) and multiplying 
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the equation (1401) on the left by A p), we find 

E^t) = z j2 [d%f^^h j (k,t)a^(0)e^-H.c. 



A=l 
2 



+ iJ2 d 3 k 



A=l 

3 



+ I E / d 3 q 



2(2tt) 3 

2(2tt) 3 
d 3 k 



"fcA 



kX 



(27T) 



C ij (uJ^k,t)K(q,k,0)e lkr '-H.c 



+ E/A 



i/=i 



d 3 A; 



(27T) 



^■(a;^ k, t)d u (q, k, 0)e lk r + H.c. 



kv 



(43) 



where = c\k\, u$ = c\q\ and we have used A(—k,p) = A*(k,p). The 
tensors 7, ( and 77 are given by 



— * 

jij(k,t) = L~ 



A(k,p) 



-1 



£ij(k, t) 



L- 1 



i(i+3) 
-1 



A(M 



P 



E [Aftp) 



P 



_p + iuj$j: 



(44) 



where L 1 {f(p)} denotes the inverse Laplace transform of f(p) and the 
tensors /, g are given by 



-ik-f 



-ik-f 
(45) 
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— * 

Similarly for the magnetic field H we obtain 



Hi{r,t) 



2 



2(2tt) 3 



3 

+ *E / dS( l 



3 



d 3 k 



kX 



kX 



(27T) 



(ij(u)$,k,t)K(q,k,0)e 



ik-r 



H.c. 



kv 





d 3 fc r 


/A/ 


,/(27r)3 ^ 



77ij(u; f , fc, k, 0)e tk r + H.c. 



kv 1 



(46) 



where 



Jij{k,t) = L 1 
lj(k,t) = -L- 1 

f)ij((jjf,k,t) = -Lr l 



-l 

(i+3)(j+3) 

-1 

(»+3)j 
3 



A(fc,p) 

p £([a(M 



P 3 



P + iUJ q l=1 



£ [A(£,p) 



(i+3)i -0' 



f,.{ujff, k) 

(i+3)(i+3) 9 7 



(47) 



The expressions f l43|) - fj47|) and the commutation relations (j3J) and (jTj) show 
that the commutation relations between components of the operators and 
if are independent of the various choices of the coupling tensors / and g 
satisfying the relations (1281) and (|29|) for a given pair of the tensors \ e an d 
X m - Finally having the fields E and H we can obtain the polarization fields 
P and M from the constitutive equations fl26l) and (T5j 



It should be pointed out that hereunto we have assumed the medium cov- 
ered by this quantization scheme is a non-conductor polarizable and magne- 
tizable medium. In the case of a non-conductor magnetodielectric medium 
the operator P defined by (jSj) is the electric polarization density and 
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as the current density induced in the medium due to the electric polariza- 
tion. When we are concerned with a conductor magnetodielectric medium 
the fields P and D in (TlBl may not be interpreted as the electric polarization 
and displacement field and the quantity ^ in Maxwell equation f|T9|) . where 
can be rewritten as 



dE dP 



V x H, 



(48) 



is not merely the current source created by the electric polarization field but 
is the sum of the free current density , due to the motion of the free charges, 
together with the current due to the polarizability of the medium. In this 
case if we compute d u (k,q,t) from fl23|) and substitute it in ([8]), instead of 
the constitutive equation ( l26l) . we obtain the linear response relation 



dPj(r,t) 
dt 

where now 
Q(r,?,t) 



J Ni (r, t) ± dt' d 3 r' Q fi (f, r', \t\ - t') E\r\ ±t'), (49) 



^ J du u 3 cos art j <f'r" 



f(u, f , r 7 ') ■ f(u , r' , r 7/ ) 



t > 

t < 
(50) 



d\ e (r,r' ,t) 



is Eq — x g t + <r(r, r', t) with x e and a respectively the electric susceptibility 
and conductivity tensors of the medium and 



Jm(r,t) 



dP, 



Ni 



dt 



3 

-*£ 

u=l 



d 3 q 



[d u (k,q,0) 



(2tt) 3 



d A k / d 6 r'u% fij(oj % ,r,r') 



(51) 



is a noise current density. In this case the explicit forms of the electromag- 
netic field clearly is the same as the relations (j4"5]) - (j4"T|) except that we must 
replace pi(k, p) by pi(k, p) + ct(k, p) where a is the transformation introduced 
in ( 1381) for the conductivity tensor a. 
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5 Summary 



By modeling an anisotropic magnetodielectric medium with two independent 
quantum fields, namely E and M quantum fields, we could investigate elec- 
tromagnetic field quantization in the presence of such a medium consistently. 
For a given pair of the electric and magnetic susceptibility tensors x e an d X™ 
of the medium, we could find the corresponding coupling tensors / and g, 
which couple electromagnetic field to E and M quantum fields respectively 
The explicit space-time dependence of the noise polarizations were obtained 
in terms of the ladder operators of the medium and the coupling tensors as 
a consequence of Heisenberg equations. In this approach, both the Maxwell 
and constitutive equations were obtained as Heisenberg equations of the total 
system. In the limiting case, i.e., when there is no medium, this approach 
tends to the usual quantization of the electromagnetic field in free space. 
Also when the medium approaches to a nonabsorptive one, the noise polar- 
izations tend to zero and this quantization scheme is reduced to the known 
quantization in such a medium. 
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